Anisotropy of tracer dispersion in rough model fractures with sheared walls. 
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Abstract 

Dispersion experiments are compared for two transparent model fractures with identical complementary rough walls 
but with a relative shear displacement ^parallel || U) or perpendicular _L U) to the flow velocity U. The structure 
of the mixing front is characterized by mapping the local normalized local transit time t(x, y) and dispersivity a(x, y). 
For S ± U, displacement fronts display large fingers: their geometry and the distribution of t(x,y)U/x are well 
reproduced by assuming parallel channels of hydraulic conductance deduced from the aperture field. For o \\ U, the 
front is flatter and a(x, y) displays a narrow distribution and a Taylor-like variation with Pe. 



Channelization is a key characteristic of flow and 
transport in fractured rocks (| 12]) and results frequently 
from the occurence of relative shear displacements of 
the two fracture surfaces during fracturation (Qi). 
Such displacements (named ^ thereafter) have been 
shown both experimentally and numerically (fill fH]) 
to create channels and ridges perpendicular to S. Their 
length depends on the multiscale geometry of the frac- 
ture walls and, even for small amplitudes S, may be a 
significant fraction of the fracture size. The permeabil- 
ity is then anisotropic: both its value and the correlation 
length of the velocity field are higher for a mean flow 
parallel to these channels (i.e. perpendicular to o). 

The objective of this communication is to demon- 
strate experimentally that this type of channelization 
also induces a strong anisotropy of the magnitude and 
properties of tracer dispersion. This is achieved by com- 
paring dispersion for mean flows parallel and perpendic- 
ular to the direction of the channels, but with identical 
flow parameters and geometry otherwise. A previous 
work ([4]) studied dispersion in one similar model (with 
a lower value of S) but with the difl'erent objective of an- 
alyzing the influence of the fluid rheology. Here, the dy- 
namics of the process, i.e. the variation with distance of 
the geometry and thickness of the mixing front is more 
specifically compared in the parallel and perpendicular 
configurations. 

Many experiments on solute spreading in fractures 
have been reported: observed dispersion co- 

efficients D increasing linearly with the mean flow ve- 
locity U (i.e. the dispersivity a = D/U is constant). 
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However, these measurements were all realized at the 
outlet of the sample with no information on the develop- 
ment of the mixing front with distance. Measurements 
by fisll used radioactive tracers, still with a resolution 
too low to investigate local spreading. In all these pa- 
pers, the anisotropy of dispersion is not investigated and 
(except for id) little information is available on the 
relative position of the fracture walls. 

We use transparent model fractures allowing for high 
resolution optical concentration measurements over 
their full area. The models are mounted vertically be- 
tween a light panel and a 16 bits Roper digital camera. 
Fluid flow takes place between two self-affine rectangu- 
lar rough walls of same characteristic exponent // = 0.8 
as in many fractured rocks (| 16]). The mean flow ve- 
locity U is parallel to the length Lx = 350 mm of the 
walls (their width is Ly = 90 mm). The two walls 
are complementary and identical in the two models and 
they match perfectly when put in contact; then, they are 
pulled away normal to their mean surface and a lateral 
shear S, parallel or perpendicular to U (i.e. to x) is in- 
troduced. In these two configurations, referred to as 
S \\ U and S ± U, the mean velocity U is therefore 
respectively perpendicular and parallel to the channels 
created by the shear. Both S and the mean aperture a are 
equal to 0.75 mm. The standard deviation of the aper- 
ture cr a = 0.144 mm is larger than for the similar models 
of fl] ((Ta = 0. 1 1mm): as a result, the flow field is found 
to be more strongly channelized. 

The fluids are shear thinning 1000 ppm solutions 
of scleroglucan in water with a high constant viscos- 
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ity (jio ~ 4500 mPa.s) at low shear rates y < 70 pre- 
venting the appearance of unwanted buoyancy driven 
flows (1 17]). One has 70 = 0.026 s"^: for a viscous 
Newtonian flow between paraUel plates at a distance 
a, the corresponding mean velocity is f/o = ^270/6 = 
3 X 10"^ mm/s. At shear rates y > 70, the viscosity de- 
creases as yu oc 7""^ with n = 0.26 (see [4]). One of the 
fluids contains 0.2 g/l of blue dye and the densities are 
matched by adding NaCl to the other. The flow velocity 
U is constant during each experiment with: 0.0024 < 
U < 0.24 mm/s and tracer transport is characterized by 
the dimensionless Peclet number Pe = Ua/Dm where 
Dm = 6.5 10"^^m^/s is the molecular difl'usion coeffi- 
cient of the dye. The experimental procedure and the 
determination of dye concentration maps from images 
recorded at constant time intervals are described by [ 4|] . 

Fig. [T] shows maps obtained in the two model frac- 
tures at two different velocities U. lfS± U, two fingers 
soar upwards with a large trough in between (Figs. [TJ) 
and[T]l): they correspond to faster paths parallel to U 
(Fig. and their amplitude increases with U. For 
S \\ U, the front is smoother (Figs.[T^ and[T]:) while its 
mean slope and the size of the indentations still increase 
with the velocity. 
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Figure 1 : Maps of the relative concentration c of the displaced fluid 
(white c = 1, black c = 0) in two transparent models, (a)-(c): S \\ U; 

(b) -(d): S ±U. Mean velocities: (a)-(b): U = 0.0125 mm/s, Pe = 14; 

(c) -(d): U = 0.25 mm/s, Pe = 285. The injected volume of displacing 
fluid is half the void space. 

The large structures in Figs. [TJ) and [TJi reffect ve- 



locity contrasts between the channels created by the 
shear. They are well reproduced by modelling the frac- 
ture aperture field as a set of independent parallel chan- 
nels of aperture a(y) =< a(x,y) >x (|1, 2]). A particle 
starting at a transverse distance y at the inlet is assumed 
to move at a velocity proportional to a(j)^"^^^/" where 
n = 0.26 for U > Uq Sindn = I for U < Uq. 

The profile Xf(y, t) of the front at a time t is then: 



Xfiy, t) 



x(t) a(3;)(»+^)/» 
< a(y)^' 



n+l)/n 



(1) 



where x(t) =< Xfiyj) >y and < aiy)^''^^''''' 
erages calculated over y. The profiles computed using 
Eq. ([T]) and the actual aperture fields appear in Figs.[Tt-d 
as dotted lines (from the above discussion, one assumes 
that n = I at the lowest velocity (Figs.[Tt and[TJ)) and 
n = 0.26 at the highest one (Figs. [It and[I]l)). Eq. ^ 
predicts well the location and shape of the "fingers" and 
"troughs" at both velocities forS±U although their 
amplitude is slightly underestimated in Fig. [T]3. In this 
latter case, one has U - Uq, corresponding to a transi- 
tion regime between the power law and Newtonian rhe- 
ologies. 

These results demonstrate that, forS±U (i.e. if U is 
parallel to the channels created by the shear), the large 
scale features of solute transport are determined by the 
contrasts between the velocities in these channels which 
increase with their aperture. Then, front spreading is 
purely convective and the total width Ax of the front 
parallel to U (i.e. the distance between the tips of the 
fingers and the bottom of the troughs) increases linearly 
with distance sls x AU/U (AU/U = typical relative ve- 
locity contrast between the different channels). These 
curves also demonstrate that the difference between the 
sizes of the fingers in the two cases are accounted for 
by the different rheological behavior of the fluids: the 
velocity contrasts (and therefore the size) are amplified 
for Pe = 285 (shear-thinning power law domain) com- 
pared to the vicinity of the Newtonian constant viscosity 
regime (Pe =14). 

For S \\ U, the features of the front are also visible at 
the same transverse distances y in Figs.[Tt and[T]:: they 
reflect again a convective spreading of the front due to 
velocity contrasts between the flow paths. However, in 
contrast with the previous case S ± U, these features 
(except for the small global slope of the front) are not 
reproduced by the theoretical model (dotted line): this 
was to be expected since its underlying hypothesis are 
not satisfied for ^ || U. The time variations of the lo- 
cal concentration c(x, y, t) on individual pixels comple- 
ment the above study by providing quantitative informa- 
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values and determining finally the distribution of the re- 
sults. As expected, the distributions are much broader 
for ^ ± ^ (Figs.[2]3-d) than for ^ || U (Figs.[2h-c). For 
Pe = 285, the experimental distribution for flow parallel 
to the channels coincides very well with the theoretical 
one and displays two peaks reflecting the structuration 
of the flow. At Pe = 14, the width and global shape of 
the experimental and theoretical distributions are overall 
similar and narrower than for Pe = 285 due to the lower 
velocity contrast in the Newtonian limit. For S \\ U 
(Figs. [2^-c), the distribution at both Peclet numbers is 
much narrower than for ^ _L U. The mean peak cor- 
responds to t(x,y)U/x ^ 1; its width increases with 
Pe (again likely due to an increase of the velocity con- 
trasts) and is similar to that of the theoretical distribu- 
tions. The experimental distribution displays however 
additional "aisles": these reflect likely complex paths 
deviating from straight trajectories parallel to U. 



Figure 2: Histograms of the experimental normalized local transit 
time t(x,y)U/x (continuous lines) for the same models and Pe values 
as in Figs.n(the letters corresponding to the different experiments are 
the same in both graphs). The distribution of t(x, y)U/x has been com- 
puted in the upper fifth of the length of the model to make meaningful 
comparisons with the distribution of the theoretical transit times ((A) 
symbols and dotted lines) determined from Eq. ([fl. 



tion on the interplay of convective and difl'usive mech- 
anisms. For all the experiments, c(x,yj) is found to 
be well fitted by solutions of a classical ID convection- 
dispersion equation for a step-like initial variation of c 
at the inlet (x = 0): 



c(x,yJ)=Ul.erfI^M^X 
2 ^/4D(x,y)t 



(2) 



Here, t(x, y) and D(x, y) are the mean transit time and 
apparent dispersion coeflftcient deduced from the time 
variation of the relative concentration c(x,yj) of the 
displaced fluid at point (x,y) (transverse dispersion is 
neglected here). It will be shown below that the two pa- 
rameters provide complementary information: D(x,y) 
(or rather the dispersivity a(x,y) = D(x,y)/U) char- 
acterizes the local thickness of the mixing front while 
t(x, y) is related to its global geometry. 

For each experiment, t(x,y) and D(x,y) are deter- 
mined for all pixels inside the field of view. Fig. [2l com- 
pares experimental and theoretical probability distribu- 
tions of the normalized times t(x,y)U/x for the same 
experiments as Fig. [T] The theoretical distribution is 
obtained by taking t(x,y)U/x = Xf(yj)/x(t), comput- 
ing the ratio Xf(y, t)/x(t) by means of Eq. ([T]) for all y 
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Figure 3: (a)-(b): Histograms (grey levels) of the values of the nor- 
malized local dispersivity a(x,y)/a (vertical scale) as a function of 
the distance x (horizontal scale) for Pe = 285. White: maximum 
probability, black : zero probability, (c), (d): Histograms obtained at 
X = 240 mm for Pe = 285 (continuous line) and Pe = 14 (dashed 
line), (a), (c): ^|| 0; (b), (d): ^±U. 

While, from the above results, the overall geometry 
of the mixing front seems to be determined mainly by 



convective effects, we examine now the relative influ- 
ence of convection and difl'usion on the local width of 
this front: this may be characterized by the variation of 
the local dispersivity a(x, y) with Pe. Figs.[3^-b display, 
for each value of x (horizontal scale), the histogram 
(coded in grey levels) of the corresponding values of 
a(x,y)/a (vertical scale). 

For ^11 U, the probability distribution of a(x,y)/a 
is narrow, particularly at high flow velocities (Figs. [3^ 
and c). Moreover, the mean value varies little with the 
distance x and reaches a constant value a(Pe)/a for x > 
100mm (see inset of Fig. (4]). The increase with time of 
the local front thickness is therefore difl'usive and can 
be characterized by a single dispersivity value a(Pe)/a. 

For S ± U (Figs.[3t) and d), the distribution of the val- 
ues of a/ a is much broader and displays a "tail" at large 
values of a /a. At the two highest velocities, the distri- 
bution displays two peaks (solid curve in Fig.O. While 
the value of a(x)/a corresponding to the peak(s) seems 
to reach a limit at long distances x, the global width of 
the distribution keeps increasing with x. In contrast to 
the case o \\ U, the increase of the local thickness of the 
mixing front with x is not simply diff'usive and cannot 
be characterized by a single dispersivity parameter. 
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Figure 4: Variation of the normalized dispersivity a/aasa function of 
Pe. ( A) ^ II U: mean value of a /a. (1)5 ± U: value corresponding to 
the first peak in the distribution of a /a. Continuous and dashed lines: 
Taylor dispersion between parallel plates respectively for a power law 
fluid of exponent n = 0.26 and a Newtonian fluid (n = 1). Inset: 
variation of the mean value < a(x, y) >y as a function of the distance 
X for 6 II U; Pe = 285 (V), Pe = 142 (A), Pe = 28.5 (□), Pe = 14 (o). 

The variation of a(Pe) with Pe for o \\ U displayed 
in Fig. 131 provides quantitative information on the lo- 



cal dispersion mechanisms: the values are similar (al- 
though slightly higher) to the predictions for Taylor dis- 
persion with a power law fluid (continuous line). At 
low Pe's, the values obtained for ^ = 1 (dashed line) 
and for ^ = 0.26 are similar and the transition towards 
a Newtonian behaviour should not influence the varia- 
tions. A dominant contribution of Taylor dispersion has 
already been demonstrated in models with a randomly 
distributed aperture of short correlation length (|7]); 
however, at low Pe's, an additional geometrical dis- 
persion regime (a = cst(Pe)) was observed and is not 
present here. 

In the opposite case (o _L U), no single value of 
a characterizes dispersion in the whole model. How- 
ever, the first peak in the distribution of a likely reflects 
dispersion in regions of low transverse gradient of the 
mean front velocity in the plane of the model: there, 
additional spreading due to transverse diff'usion in the 
velocity gradient should be reduced. The values of a/ a 
corresponding to this first peak have therefore been plot- 
ted in Fig|4](T) as a function of Pe. At high velocities, 
they are indeed close to those corresponding to Taylor 
dispersion (and to the other model); at low velocities, 
they remain higher. 

Overall, the present experiments demonstrate that so- 
lute dispersion in a channelized rough fracture depends 
crucially of the orientation of the flow. Optical mea- 
surements allowed us to characterize the dynamics of 
dispersion at both the global and local scales from mean 
transit times and local dispersivities deduced from local 
concentration variations. 

For ^ ± U (U parallel to the channelization), the 
large scale geometry of the displacement front is con- 
trolled by the velocity contrasts between the channels. 
At all Pe values (14 < Pe < 285) the geometries of 
the large fingers and troughs in the front and the dis- 
tribution of the local transit times t(x,y)U/x are well 
predicted from a transverse eff'ective permeability pro- 
file computed using the aperture field. Moreover, dif- 
ferences between the finger sizes at low and high Pe's 
are explained by variations of the fluid rheology (shear- 
thinning at high Pe's and Newtonian at low ones). These 
results confirm fully the convective origin of the large 
scale structures of the front and the relevance of local 
measurements all along the flow paths: measuring only 
variations of the mean concentration at the outlet (as is 
often done practically) and fitting it to a solution of a 
convection-diff'usion equation might indeed lead to an 
incorrect identification of such processes as geometri- 
cal dispersion. 

For II U (U perpendicular to the channelization), 
the front is flatter and the distribution of the transit times 



is narrower than for ^ _L IJ\ this reflects a more efl'ec- 
tive sampHng of the velocity heterogeneities by the so- 
lute particles. As could be expected, the remaining geo- 
metrical features of the front cannot be predicted by the 
channel model: however they are still observed at the 
same transverse locations at all velocities and remain 
of similar amplitudes (taking into account the variations 
of the rheology). This suggests that these features are 
again of convective origin. 

Still for 6 \\ U, the variation of the local front thick- 
ness is difl'usive and well characterized by a single mean 
dispersivity a(Pe); its dependence on Pe suggests a 
dominant influence of Taylor dispersion. For o ± U,m 
contrast, the broad distribution of a(x, y) does not allow 
one to define a meaningful global dispersivity. In ad- 
dition, even the lowest values of a(x, y) (corresponding 
to simple flow paths) are larger than those expected for 
Taylor dispersion except at the highest Pe values; values 
corresponding to Taylor dispersion are also expected at 
smaller values of S leading to a reduced disorder of the 
flow field f7] . This latter diff'erence, as well as the tail in 
the distribution, reflect the increasing influence of trans- 
verse molecular difl'usion inducing tracer exchange with 
adjacent flow paths of difl'erent velocities and enhancing 
dispersion. 

This set of results is highly relevant to the interpreta- 
tion of field observations (p"]). However, the influence 
of the length of the samples is an important issue and 
will need to be investigated before transposing these re- 
sults. In particular, although the front geometries are 
determined in both configurations by spatial variations 
of the flow velocities, a geometrical dispersion regime 
has never been observed, in contrast with experiments 
on rough fractures with a small correlation length of 
the aperture (itzl]). For much longer fractures, trans- 
verse difl'usion might be large enough for solute parti- 
cles to sample the whole distribution of local velocities 
and reach such a global difl'usive spreading regime. In 
the case of a broad distribution of the hydraulic conduc- 
tivities, one might also observe instead an anomalous 
dispersion regime as suggested by |5] and |6]. Because 
of the specific correlations of the flow velocity field for 
self-affine wall geometries (Jll]), the corresponding ex- 
ponent would then likely depend on the characteristic 
roughness exponent of the fracture walls. 

1. Acknowledgments 

We are indebted to R. Pidoux for his assistance. HA 
and JPH are supported by CNRS through the GdR No. 
2990) and by the EHDRA (European Hot Dry Rock 
Association). This work was greatly facilitated by a 



CNRS-Conicet Collaborative Research Grant and by 
the Ecos Sud A03E02 program. 

References 

[1] Auradou, H., G. Drazer, A. Boschan, J. P. Hulin, and J. Koplik 

(2006) , Shear displacement induced channelization in a single 
fracture, Geothermics, 35, 576-588. 

[2] Auradou, H., A. Boschan, R. Chertcoff, S. GabbanelH, J.P. 
Hulin and 1. Ippolito(2008), Enhancement of velocity contrasts 
by shear-thinning solutions flowing in a rough fracture, /. Non- 
Newtonian Fluid Mech., 153, 53-61. 

[3] Becker, M. and A. Shapiro (2000), Tracer transport in fractured 
crystalline rock: Evidence of nondifl'usive breakthrough tailing. 
Water Resour. Res., 36, 1677-1686. 

[4] Boschan A., H. Auradou, I. Ippolito, R. Chertcofl' and J.P. Hulin 

(2007) , Miscible displacement fronts of shear thinning fluids in- 
side rough fractures. Water Resour Res., 43, W03438. 

[5] Bouchaud J.P., A. Georges and P. Le Doussal (1987), Difl'usion 
anormale dans les milieux desordonnes : piegeage, correlations 
et theoremes de la limite centrale, /. Phys. France 48, 1855- 
1860. 

[6] Dentz M., Le Borgne T. and Carrera J. (2008) Efl'ective transport 
in random shear flows Phys. Rev. R 77, 020101. 

[7] Detwiler R., H. Rajaram, and R. J. Glass (2000), Solute trans- 
port in variable aperture fractures: An investigation of the rela- 
tive importance of Taylor dispersion and macrodispersion. Wa- 
ter Resour Res., 36, 1611-1625. 

[8] Gentler, S., E. Lamontagne, G. Archambault, and J. Riss (1997), 
Anisotropy of flow in a fracture undergoing shear and its rela- 
tionship to the direction of shearing and injection pressure. Int. 
J. RockMech. &Min. Sci., 34(3-4), 412. 

[9] Keller A.A., Roberts RV., Blunt M.J. (2000) Efl'ect of fracture 
aperture variations on the dispersion of contaminants. Water Re- 
sour Res., 35, 55-63. 
[10] Lee J., J.M. Kang and J. Choe (2003) Experimental analysis on 
the efl'ects of variable apertures on tracer transport. Water Re- 
sour Res., 39 (1), 1015, doi : 1 0. 1 029/200 IWROO 1246, 2003. 
[11] Matsuki K., Y Chida, K. Sakaguchi and RW.J. Glover (2006), 
Size efl'ect on aperture and permeability of a fracture as esti- 
mated in large synthetic fractures. Int. J. Rock Mech. Min. Sci., 
43, 726-755. 

[12] Committee on Fracture Characterization and Fluid Flow, Na- 
tional Research Council (1996), Rock Fractures and Fluid 
Flow: Contemporary Understanding and Applications, Na- 
tional Academy Press, Washington, D.C. 

[13] Neretnieks, I., T. Eriksen, and P. Tahtinen (1982), Tracer move- 
ment in a single fissure in granite rock: Some experimental re- 
sults and their interpretation. Water Resour Res., 18, 849-858. 

[14] Olsson, W. A., and S. R. Brown (1993), Hydromechanical re- 
sponse of a fracture undergoing compression and shear. Int. J. 
RockMech. Min. Sci. & Geomech. Abs., 30(1), 845-851. 

[15] Park C.K., T.T. Vandergraaf, D.J. Drew and RS. Hahn (1997), 
Analysis of the migration of nonsorbing tracers in a natural frac- 
ture in granite using a variable aperture fracture model, /. Cont. 
Hydro. 26, 97-108. 

[16] Poon C, R. Sayles, T. Jones, Surface measurement and fractal 
characterization of naturally fractured rocks, J. Phys. D: Appl. 
Phys., 25 (1992) 1269-1275. 

[17] Tenchine S. and Ph. Gouze (2005) Density contrast eff'ects on 
tracer dispersion in variable aperture fractures. Advances in Wa- 
ter Resources, 28, 273-289. 



